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Abstract. We study a class of superconductive radiation detectors in which the 
absorption of energy occurs in a long superconductive strip while the redout stage is 
provided by superconductive tunnel junctions positioned at the two ends of the strip. 
Such a device is capable both of imaging and energy resolution. In the established 
current scheme, well studied from the theoretical and experimental point of view, 
a fundamental ingredient is considered the presence of traps, or regions adjacent to 
the junctions made of a superconducting material of lower gap. We reconsider the 
problem by investigating the dynamics of the radiation induced excess quasiparticles in 
a simpler device, i.e. one without traps. The nonequilibrium excess quasiparticles can 
be seen to obey a diffusion equation whose coefficients are discontinuous functions of 
the position. Based on the analytical solution to this equation, we follow the dynamics 
of the quasiparticles in the device, predict the signal formation of the detector and 
discuss the potentiality offered by this configuration. 
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1. Introduction 

Superconducting single photon radiation detectors, such as transition edge sensors 
microbolometers, superconducting tunnel junctions, kinetic inductance detectors, and 
superconductive resonator detectors are attractive devices because they have nearly 
achieved the desired performance for energy-resolved detection of individual photons or 
particles over a broad energy range with high counting rates [1, 2]. 

Current work is generally directed towards two objectives: improving energy 
resolution and scaling to cover large areas. For microbolometers the current trend is to 
try to multiplex an array of detectors to one readout system. Various techniques have 
already demonstrated their viability, although several difficulties remain to be solved. 
The multiplex approach is much more difficult in the case of STJ detectors, as the signals 
are much faster, thus increasing the bandwidth necessary for the multiplexing circuit. A 
possibility to simplify the problem is to drastically reduce the number of channels that 
has to be readout without reducing the number of pixels and the area covered by the 
array. This can be accomplished by substituting several single pixel STJs in the array 
with long and narrow strips where radiation is absorbed (as long as possible to cover the 
largest area). Each strip is readout through two STJs placed at the ends[3, 4, 5]. Such 
a structure is sometime referred in literature with the acronym DROID which stands 
for Distributed Read Out Imaging Device [6]. In DROIDs energy resolution is obtained 
in the same way as for a single STJ by counting the total number of quasiparticles that 
tunnel through the barriers of the two junctions in time coincident pulses; the position 
is obtained by measuring the quasiparticles separately collected by the two junctions. 
The DROID scheme has been demonstrated to give very good position resolution along 
the absorber strip, and also a good energy resolution [7]. 

Since the very beginning when DROIDs were proposed, they were designed by 
using a higher gap superconductor for the absorber and a lower gap superconductor 
for the STJs electrodes [3]. The idea was to exploit the quasiparticle trapping and the 
multiplication effects [8] to improve the collection efficiency. In a DROID quasiparticles 
produced in the absorber after photon absorption will diffuse toward the opposite sides 
where they are soon trapped and then counted by the tunnel junctions. There are 
currently two different ways on how to realize such traps on the sides of the absorber. 
One places the trap on top of the absorber [3, 6, 9, 10], the other places the trap 
laterally to the absorber [5, 11]. The first method has the advantage of the best possible 
interface between the absorber and the trap material as they can be both deposited 
in situ without breaking vacuum during fabrication, but constrains the volume of the 
trap and junction area severely, and results in STJs affected by the proximity effect. 
The second has almost no limitations on trap design, and junction area, the STJ are 
not affected by the proximity effect, but much care must be taken in fabrication to 
keep the interface between the trap and absorber defect-free. In both designs energy 
resolution has been measured as degraded because of the temperature rise of the trap 
when a large amount of quasiparticles enter [12, 13]. This degrades the ability of the 
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trap to efficiently capture the quasiparticles and increases the dark current of the ST J, 
which results into a degraded position and energy resolution. The easiest way to remedy 
this problem is to make bigger traps, if possible, which would reduce the heating effect 
of the quasiparticles. Apart from the case where the trap volume is constrained by 
the absorber design, larger traps also means that losses and diffusion inside the traps 
become relevant phenomena that have to be considered. 

In the light of the above considerations, we found the physics of a different DROID 
design, one that does not use traps, interesting to investigate. The presence of traps 
is not a fundamental ingredient that underlies the DROID idea. The important 
issues are the large area coverage with imaging capability and the reduction of the 
readout channels, while maintaining good energy resolution. Moreover, single pixel STJ 
detectors have been made without traps; they are well proven examples of detectors that 
does not have traps but conserve excellent energy resolution [14]. Large area coverage 
and reduction of readout channels can be kept by simply connecting two STJs to the 
same bottom electrode (the long absorber) without introducing traps. 

A no-trap DROID structure was already analysed in the past [15], but it was not 
fully appreciated for the advantage of being a simpler structure from several points of 
view: fabrication technology, diffusion, self-heating. In the model contained in [15] the 
junctions were considered point-like and were only sensitive to the distribution function 
of quasiparticles at the edges of the device without any disturbance of the nonequilibrium 
state in the absorber. In the present work the analysis has been extended to consider the 
finite size of the junctions, and to include their influence on the nonequilibrium dynamics 
of excess quasiparticles. In this way the presence of the junctions can determine the 
detector response. 

Finally we remark that the investigation of a DROID structure without traps is also 
of interest in the perspective, just recently emerged[16], that to obtain STJ-detectors 
which would be really competitive in terms of energy resolution one will be probably 
forced to use for the material of the absorber a superconductor with energy gap equal 
or lower than that of aluminium (A < 0.18 meV). The reason is that the ultimate 
energy resolution in STJs is proportional to the square root of the energy gap of the 
superconductor where radiation is absorbed. Keeping an energy E = 1 keV as reference 
for X-ray photons, the intrinsic energy resolution requirement for advanced applications 
is 2 eV. This is achieved when the energy gap value is equal or less that of the aluminium. 
Moreover the present state-of-the-art in STJ fabrication technology is based on the 
A1 2 3 oxide tunnel barrier, indicating again STJs with Al electrodes as the best and 
possible choice. It is then natural to consider a whole aluminium DROID, an Al-absorber 
laterally readout by Al-STJs: definitively a no-trap DROID structure. 

In this paper we describe the response of such a detector under the influence of 
single photons. A model is developed, and an analytic expression is obtained for the 
evolution in time of the one-dimensional quasiparticle density, from which we calculate 
the tunnel currents and the collected charge through both the STJs. We discuss the 
position dependence of the pulse height governed by diffusion and quasiparticles decay. 
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In a superconductor radiation energy is absorbed in a complicated process that results 
in the generation of an excess number of quasiparticles [17]. A DROID measures the 
number of nonequilibrium excess quasiparticles as they tunnel through the tunnel junc- 
tions positioned at the end of the superconductive absorber. Each junction is polarised 
at a convenient voltage point in the subgap region of the I-V characteristic. When 
illuminated, the device provides energy and position information of the incoming radi- 
ation by comparing the coincident measured charges in the two junctions. To describe 
the dynamical behaviour of the proposed detector a model of the spatial and temporal 
quasiparticle excess density is needed. The structure is schematised in figure 1. It con- 
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Figure 1. Schematic of the cross section of the proposed detector (not in scale). 
Dashed arrows on the two sides indicate the excess quasiparticle tunnelling 

sists of an absorber of length L, with tunnel junctions at the two edges whose length 
is (L — a)/2. We assume that the effect of the tunnel junctions is to remove the excess 
quasiparticles created in an impact event of a photon with the absorber with a constant 
rate probability, 7t U nn, and that the quasiparticles do not re-enter the absorber once 
removed. We also assume, in the absorber, isotropic diffusion constant, D, and uni- 
form constant loss rate, 7i OS s, which accounts for the various mechanism of quasiparticle 
losses. These assumptions enable us to describe the spatial and temporal evolution of 
the excess quasiparticle density, n(x,t), through the standard diffusion equation [4] 

dn(x,t) B 2 n(xA) , . . . T . 

m = D d 2 -l(x)n(x,t), -L/2 < x < L/2, < t < oo (1) 

with the boundary conditions 

dn(-L/2,t) = dn(L/2,t) = Q 
dx dx 

The presence of the tunnel junctions is then modelled through 7 (x) = 7, where the 
index i = 1, 2, 3 refers to the three regions defined as 
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{7i = 7ioss + Ttunn if -L/2 < x < -a/2 Region 1 
72 = 7ioss if -a/2 < x < a/2 Region 2 (3) 

73 = 7io SS + 7tunn if a/2 < x < L/2 Region 3 

Thus in the two junction regions the probability of tunnelling adds to the probability 

of loss. Although simple looking, this problem is more complicated compared with 

that arising when traps are present and, as far as we know, the solution has never been 

reported. We choose to solve this problem by splitting it into three parts, corresponding 

to the regions 1, 2, and 3 above. Thus we need to supplement the boundary conditions, 

Eq. 2, with the requirements that the excess quasiparticle density, and its spatial 

derivative, are continuous functions at x — ±a/2 

f a \ ( a \ (a \ i 'a \ 
n 1 (--,t)=n 2 (--,t),n 2 (-,t)=n 3 (-,t) 

dm (-It) _ dn 2 (-f,t) <9n 2 (f,t) _ <9n 3 (f,t) (4) 
dx dx dx dx 

The radiation impact on the absorber is modeled through the initial condition: 

n (x, 0) = N 5 (x - s ) , (5) 

where x is the radiation impact point and N the number of quasiparticles generated 
by the radiation impact. All this provides seven relations which together with equation 
(1) allows to determine the solution in the three domains between —L/2 and L/2. 



3. Solution of the problem 



The solution to the problem, equations (l)-(5), can be obtained by the separation of 
variables method [18] introducing the arbitrary constant v. Then the unknown density 
can be written as an integral on all possible values of v in the three previously defined 
spatial domains. 



rii (x, t) 



Xi (ui) e 



+ 



-(7»+"»)t 



Xi (vi) cos \l jjX + Yi (u^ sin 



-x 



dvi. (6) 



Inserting the full formal solution, equation (6), into the six boundary conditions, 
see Appendix A, we find that the possible values of v is an infinite discrete set. It is 
convenient to distinguish the sign of v by introducing a 2 = —Dv/L 2 and fi 2 = Dv/L 2 . 
Then the two equations determining these values of v, written in terms of a and fi 
respectively, are 



[7d + (td — 2a 2 ) cosh (a — aa)~\ sin ^oa/7d — « 2 j + 
— 2q;a/7d — o 2 cos ^<2a/7d — o 2 ^ sinh (a — aa) = 



(7) 
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and 

[7d + (td + 2yu 2 ) cos (// - a/i)] sin (a^/j D + /x 2 ) + 

+2/iA/7o + /i 2 cos (a\/7D + At 2 j sin (// - a/i) = 0, (8) 

where 7d = 7t un n-^ 2 /-D- Once the possible values of the constant v have been determined 
the two integrals in equation (6) change into sums. Thus equation (6) can be written as 

(x,r) = Rn(oc n , x )e UoL Q "^ r F iin (a n , x) + 



N a —1 

Ui(x,T) = 

n=0 

+ J2 S n (fi n ,x )e-(^ + ^) T G hn (fi n ,x) (9) 

n=JV a 

where a n , fi n are solutions to equations (7) and (8), respectively, and the two 
proportionality constants R n (a n ,xo) and S n (fi n ,xo) will be determined by the initial 
condition. We have also introduced the normalised time, r = t7i OS s, the finite number 
of possible a- values, N a (see Appendix A), the diffusion length A = a/-D/7i OS s, and the 
spatial base functions 

f An K) e°" x + fl i>n (a n ) e~ a " x , i = 1, 3 

X [ A,n («n) COS (V7D - "n 1 ) + 5 i,n («n) sin (a/V) ~ « 2 2;) , % = 2 

f Ci,n COS (/i n x) + A,n (/in) Sin (/i n x) ^ = 1,3 

' " / '""' I C i>n (/!„) cos (a/7d + /i 2 a;j + D i>n (// n ) sin (^a/td+T^^) , i = 2 

in which all lengths (x,x ,a) are normalised to L, unless otherwise indicated. The 
procedure to obtain the coefficients A iin , B itn , C^ n , D itn , R n , S n is described in detail in 
Appendix B along with their analytical expressions. 

At this point the mathematical problem is completely solved, and the spatial 
and temporal response of the proposed detector is known. An example of the spatial 
evolution of the density of excess quasiparticles in an aluminium device is shown in figure 
2 at subsequent instants of time (in units of 7i OS s) for an event occurred at x = 0.3. 
We used values of the parameters which are compatible with those of a device based on 
aluminium: D = 60 cm 2 /s, 7i oss = 10 4 s -1 , A = 775 /im, L = 775 /im, for the absorber, 
and a tunnelling rate 7t U nn = 1-43 x 10 6 s^ 1 for the two junctions, each long 50 /im, 
(a = 675 /im, in absolute units). Note that no trace remains, in the time evolving 
density, of the discontinuity in the coefficient 7(2;) so that the transition between the 
absorber and the junction regions is completely smooth. 

4. Analysis of the detector signal 

To show the potential of the result obtained in the previous section, we will calculate 
the signals that are typically measured in experiments with DROIDs (i.e. tunnelling 

current pulses flowing through the junctions). These currents are calculated as 

„_« 

h (x , t) = 7tun n e-L / n (x, r) dx (12) 
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Figure 2. Temporal evolution of the density of excess quasiparticles in the device 
for an event occurring at x = 0.3. Curves 1,2,3,4 correspond to r = t7i OS s = 
0.0005, 0.001, 0.01, 0.5. The parameters used are A/L = 1, 7 tU nn/7io SS = 143, a = 0.87 



h (x , t) = 7t mm eL / n (x, r) dx 
and are given by the expressions 



(13) 



t / \ r ^V^ 1 Rn (OL n , Xo) -(-Oil a 2±i) T { a n acc n \ / l+o \ 

h (x , t) = 7t u „ n eL ^ V 7 e U«. "^J fe" - e~ ) U ljn e a " — + B ljn ) 

„ =0 «n 



n=N a 



X 



COS 



a/i r , 



sin 



sin 



} (14) 



x {fl 1>B [cos 

/ 2 (ar , r) = 7tunn eL N f t^^A^-^y (e* - e-) (^e- 3 ** + 



?1=0 



tmm eL jr s ^ Xo K-(^ + ^> 

n=N a /in 



X 



3,n 



[ COS I" 



cos 



/J, 



+ C3,ra 



sin 



/'■I, 



sin 



(15) 



Figure 3 shows the two current pulses for the absorption event shown in figure 2. In 
literature one typically finds expressions for the collected charges at the two lateral sides 
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Figure 3. Temporal evolution of the normalised currents passing through the tunnel 
junctions for the absorption event located at x — 0.3. Parameters correspond to 
A/L = 1, 7tu„n/7ioss = 143, a = 0.87 



of the absorber, which were directly derived through equations governing the charge 
itself [4] rather than the density. We emphasise that to fully characterise the DROID 
behaviour it is certainly of great help to also have expressions for the extra tunnelling 
currents as those calculated above, which is possible only when the solution for n(x, t) 
has been obtained as we did here. 

Finally the collected charge in each tunnel junction as a function of the impact 
point x can be calculated as 



Ql,2 (Xq) 



and is given by the expressions: 

N a -1 



Q\ M = 



7tu 



h,2 (x ,t) dt. 
s: 

Rn (a n , Xq) 



7loss 



=0 Ol r , 



71 

7lose 



2 A 2 
<17 



e 2 



"Qn \ / . „ l+o _ \ 

— J [A^ n e- a -— + B hn j 



7t« 



7ic 



S n (fJ'n, Xq) 



n=N a l^n 



71 



M'n L 2 J 



{ D l,n C0S (y 



COS 



apL r 



+ 



Q2 (x ) 



7tu 



N a -1 



7loss 



R n (a n , Xq) 



( an o«n \ /_ __, 1 + n . \ 



(16) 
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. 7tunn T S n (/i n ,x ) r r fa/i n \ /AM" 
H eL > t r- < D-i n cos — cos — + 

+C3,„[ S in(^)- S in(^)]} (18) 

Through equations (17), (18) curves at constant energy and constant absorption positions 
can be constructed (see next section) for a direct comparison with the typical 
experimental data representation [3, 4]. 

5. Results and discussion 

The behaviour of a DROID is often analysed by plotting coincident collected charges 
(Qi, Q-i) for all the possible values of xq at constant energy (i.e. constant No). Following 
this choice, in figure 4a we show the response of the no-trap-DROID (equations (17) 
and (18)) for various values of A/L corresponding to absorbers with different values 
of the diffusion constant D or different values of the loss rate 7i oss - The initial and 
the final points of each curve, where charge maximises, correspond always to the —L/2 
and L/2 positions on the absorber. Therefore spatial resolution scales inversly with 
charge collected in these points. It is seen that increasing A/L, simultaneously increases 
the charge yield, makes the charge yield more uniform (the curve shape becomes more 
straight) and reduces the spatial resolution (the total extension of the curve shortens). 
In the limit of very efficient diffusion (curve 4), the spatial resolution is lost and the 
response becomes similar to that of a single STJ. Figure 4b shows the no-trap-DROID 
response for increasing values of 7t U nn/7ioss> corresponding to devices with increasing 
tunnel junction transparencies. In this case the charge yield, the uniformity of the 
charge yield and spatial resolution of the device simultaneously increase. Curve 3 in 
figure 4a and curve 4 in figure 4b coincide and refer to the same device considered in 
figure 2 and 3. We stress that these curves predict a behaviour comparable to that of 
a DROID based on the quasiparticle trapping principle [4, 5] and that all the essential 
features of this last device are conserved in the device described here. 

It is also of interest to briefly discuss the role of the junction size by investigating the 
role of the geometrical parameter a. This can be done by using the same type of Q\, Q2 
plot as before by simply emphasising points corresponding to the same equidistant 
positions of radiation absorption. Two such curves are shown in figure 5 in which 
the open circles, superimposed to the solid line, mark the points corresponding to an 
arbitrary segmentation of the device in twenty equal intervals. The lower curve, curve 
1, has a = 0.95, the higher, curve 2, has a = 0.5, moreover 7tunn/7ioss = 143, A/L = 1; 
the geometry of the devices corresponding to these curves are sketched in the inset of 
the same figure (the gray regions indicate the junctions). A smaller a value (larger 
junctions) results in an increased charge yield with larger uniformity. Furthermore the 
spatial resolution is strongly nonuniform as indicated by the lateral crowding, and the 
simultaneous central rarefaction of points. Therefore the DROID becomes very position 
sensitive in the central region while under the junctions this sensitivity is degraded. To 
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normalized charge Qj/CNq e) 



1 




normalized charge Qj/CNq e) 



Figure 4. a) Collected charges Q2 vs Qi for increasing values of 
A/L, (0.25,0.5, 1.0,5.0), 7tunn/7ioss = 143 and a = 0.871.6) Collected charges Q 2 vs Qi 
for increasing values of 7 tU nn/7io SS , (10.0,25.0,50.0,143.0,400.0), A/L = l,a = 0.871 
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better illustrate this effects, in figure 5, we indicated two points (solid circle) on each 
curve corresponding to the same impinging radiation position (arbitrarily chosen) on 
the two detectors and marked in the inset by arrows. 

In conclusion we have analysed in detail the behaviour of a class of superconductive 
radiation detectors, denominated DROID, in which the absorption of energy occurs in 
a long superconductive strip and the readout stage is provided by superconducting 
tunnel junctions positioned at the two ends of the strip. We have introduced boundary 
conditions suitable for the absence of trapping regions and solved analytically the 
resulting diffusion problem. We underline that the model provides not only the total 
charges collected by the two junctions but also explicit expressions for the tunnelling 
current pulses. Our analysis shows that a device based on a single superconductive 
material (i.e. without traps) is still capable both of position and energy resolution. 
Recently, aluminium has been proposed as the material of choice to achieve the best 
energy resolution and imaging is obtained by arrays of single superconducting tunnel 
junctions. Unless other considerations are put forward, our work demonstrates the 
feasibility of an alternative viable solution to imaging, i.e the DROID configuration, 
which can be entirely based on a single material, e.g. aluminium. 




normalized charge Qj/CNq e) 



Figure 5. Detector response for two different values of the geometrical parameter a 
(the size of the junction-free region of the device). Curve 1, a device with a = 0.95, 
curve 2 a device with a = 0.5. In the inset, a sketch of the two detectors corresponding 
to the two curves. The gray regions indicate the junctions. Open circles indicate 
signals coming from equi-spaced and uniformly distributed positions along the device. 
The two solid circles are signals coming from the position marked by the arrows on 
the two devices 
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Appendix A. Matrix coefficients 

Inserting the full formal solution, equation (6), into the six boundary conditions, 
equation (2), and using the substitutions V\ — v 3 — v and u 2 = v + 7 tU nn, it is 
found that time-independence of the boundary conditions splits the set of equations 
for Xi, Y"i, X 2 , Y 2 , X 3 , Y 3 into three sets, depending on the value of v. The first set of 
equations, for the values — 7 tU nn — 7ioss < v < — 7tunn, arises from matching exponential 
functions in region 1 and 3 with exponential functions in region 2. The second set of 
equations, for the values — 7 tU nn < v < 0, arises from matching exponential functions in 
region 1 and 3 with the trigonometric functions in region 2. The last set of equations, 
for positive values of u, arises from matching trigonometric functions in region 1 and 3 
with trigonometric functions in region 2. Despite the differences between the three sets 
of equations it is possible to write them in the same matrix form 



mn 


m 12 
















" x 1 ' 




' " 


m 2 i 


m 2 2 


m 23 


m 24 

















m 31 


m 32 


m 33 


m 34 








X 


x 2 













m 43 


777.44 


77745 


777 46 


Y 2 













m 53 


777 54 




777 56 




x 3 



















r77 65 


77766 




Y 3 




_ 



(A.l) 



The exact expression of the matrix elements m^j depends on the three possible ranges 
of values of v. They are listed in Table (Al), written in terms of a 2 = —Du/L 2 and 
/i 2 = Dv/L 2 . Inspection of the determinants of the three matrices shows that for 
— 7tunn - 7ioss < v < — 7tunn only the trivial solution exists. On the contrary, non- 
trivial solutions can be found for a finite number, A^ a , of values of v when the interval 
— 7tunn < v < is considered (solutions to equation (7)), and also for an infinite discrete 
set of values of v for v > (solutions to equation (8)). These two equations express the 
conditions for theexitancee of non-trivial solutions in the two considered intervals of v. 
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Table Al. Expressions in terms of a and \i of the matrix elements 
appearing in cq. (A.l). (a 2 = -Dv/L 2 and fi 2 = DvjL 2 ). 



Element —7 — j t < v < —j t —-y t < v < 



v > 



m 12 
m 21 

"7 22 

m 23 

m 24 

m 31 
m 32 
m 33 
m u 

77743 
777.44 
77745 

77746 

m 53 

777 54 
"7 55 

m 56 

"765 



exp (-f ) 
-exp (f) 
exp (-f ) 
exp (f ) 



-7c 



exp 



— exp 

a 77721 
-a77722 

y/ct 2 ~ 77J777 23 
- a/« 2 ~ 7 £>"7 24 
a/o 2 - 77J77724 



-ID 



^ - 71)^23 
a: 77722 
— Q;77721 
77724 
"7 23 
"7 2 2 
77721 
-77712 



exp (-f ) 

-exp(f) 
exp (-f ) 
exp (f) 

ay/- 



— COS 



sin 



iD-a 



am 2 i 
-am 22 

-\flD 



a 2 777 2 4 



y /lp 

\/lD- 



a 2 m 23 
a 2 m 24: 



a/TO - a 2 777 2 3 

a77722 

— a777 2 l 

"723 

-77724 

"7 2 2 

"7 2 1 

-77712 

-mn 



sin (f ) 

cos (f ) 

COS (f ) 

-sin(f) 

' COS 2 




-/U77722 

/7"721 

~ a/7D + / 7 2 "7 2 4 
\/lD + A* 2 7"23 



/i77722 

/777721 

"723 

-77724 

"721 

-77722 

-mn 

m 12 



H 2 m 2i 



(i 2 m 23 
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Appendix B. Coefficients in the solution, equation (9) 

As the problem is symmetric in space the spatial modes separate in even and odd modes. 
The first mode is even, the second odd and so on. This means that the symmetry of 
the first mode in the second sum in equation (9) depends on whether N a is even or 
odd (which is why we have chosen to let is start from N a instead of from 0). Insight 
into when a mode leaves the second sum in equation (9) and appears in the first sum 
is obtained by examining equation 8 in the limit fi — > 0. It is found that equation 8 
reduces to sin (a^/Trj) = 0. From this limit it is easily seen that it is possible to express 
N a as the smallest integer greater than a^/ju/'K. 

We obtain exact expressions for the even coefficients (A i>2n and B i2n ) by assuming 
A 2t 2n known. Then we create a reduced system of five equations with five unknowns 
by moving all the terms containing A 2)TL to the right hand side in equation (A.l) 
and eliminate a row (exactly which one is not important as the six rows arelinearlyy 
dependent on each other). By solving this system of equations we obtain A i)2n and 
B i2b except for a proportional constant, R n (a n ,x ) in equation (9), which will be 
determined using the initial condition (equation (5)). Care should be taken as the 
expressions obtained are not valid for the odd coefficients (Ai >2n+ i and Bi j2n+ i) because 
for these coefficients A 2t2n+ i is zero and thus creates problems when selected as known. 
We instead select B 2>2n+ i as known and repeat the procedure as above. Likewise for the 
second sum in equation (9) we select first C 2j2n as known to obtain the even coefficients 
and then select D 2:2n+ i as known to obtain the odd coefficients. Also for the these 
coefficients we have to determine a proportional constant, S n (fi n ,x ) in equation (9), 
using equation (5). The expressions obtained are: 




(B.l) 



(B.2) 



(B.3) 



(B.4) 



A 




[exp (a n ) + exp (aa n )\ y/j D - a\ cos 



( 




D 



2 



QL 




+ (B.5) 



+ [exp (a n ) - exp (aa n )] a n sin 




B 2 , n (a n ) 



[exp (a n ) — exp (aa n )] a n cos 



- [exp (a n ) + exp (aa n )\ y/^D- «n sin 




(B.6) 
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C 2 ,n (/in) = V ID + K COS I I COS I — ^ I + (B.7) 



y/lD + /4 



H n sm | ] sin 



D 2 ,n (/in) = - /in Sin ( 1 COS I — V - - ) + (B.8) 



/ - //i n -a// n \ . / ay/-f D + nl 

V Id + Hi cos sm 1 



A 3 , n (a n ) = (-1)" exp (^) (B.9) 

B 3 , n («„) = (-l) n exp (a n + ^) \JlD~OL n (B.10) 

^3,n (/in) = ("I)" V^T74 COS (^) (B.ll) 

£>3,n (/in) = (-1)" v^+^sin (y) • (B.12) 

Note that if the second sum in equation 9 is started at zero instead of N a then the 
coefficients C 3 ^ n and D 3n should be multiplied with {—l) Na . 

The initial condition, equation (5), can be used to find R n (a n ,x ) and S n (fi n ,x ) 



as: 



D / \ at Fi n \ a ni x o) i \ i\t ^i,n (/^nj Xq) . . 

Rn \OL n , X ) = N —— -2—, S n {fin, X ) = Nq—— -5— (B.13) 

1 a {Ctn) f n (/in) 

Where the norm of the modes, 

gO n (l + a) 

T Q (a n ) = L [(2 - a) 7 D a n - 2a„ 3 + a7 D a„ cosh (a n - a n a) + 

+27 D sinh (a n - a n a)\ (B.14) 

and 

(/in) 2 = ~~7~~ [(2 - a) 7D/in + 2/i„ 3 + a7 D /i n cos (//„ - a/i n ) + 27^ sin - a/i„)] (B.15) 

4/in 

have been calculated using the equations (7) and (8) in order to simplify as much as 
possible their form. To obtain equation (B.13) we have used the orthogonality of the 
modes when integrated over the entire length. The exact expression of equation (B.13) 
depends on whether x is located in region 1, 2, or 3. 



